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The  Llixha  Elbasan  hot  springs  are  among  the  most  important  geothermal  springs  of  Albania.  Balneological 
use  dates  back  centuries,  but  the  first  modern  use  started  in  1937.  Unfortunately  this  water  has  not  been 
used  for  its  energetic  values  yet.  The  temperature  of  the  water  is  above  60  °C  and  the  flow  above  161/s, 
thus  direct  utilization  is  possible,  in  particular  for  space  heating.  Three-dimensional  temperature  field 
calculations  and  engineering  calculations  on  a  heating  system  with  heat  exchangers  are  presented  here. 
The  results  show  that  the  water  temperature  is  expected  to  be  stable  and  considerably  higher  temperature 
is  expected  through  deep  well  drilling.  Geothermal  water  of  Llixha  hot  springs  fulfils  all  requirements 
needed  for  a  district  heating  system  in  the  region.  The  University’s  Campus  of  Tirana  is  composed  of  29 
buildings,  which  are  partially  heated  through  a  coal  heater.  The  installed  capacity  is  2558  kW  while  the 
coal  consumption  is  about  920  kg/h.  The  University’s  Campus  of  Tirana  is  one  of  the  most  important  areas 
and  with  the  highest  density  of  population  in  Tirana,  so  it  is  the  best  area  to  show  the  heat  exchanger 
efficiency.  The  economical  analyses  prove  that  the  borehole  heat  exchangers  are  more  convenient  than 
the  coal  heating  systems. 
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1.  Introduction 

Albania  is  a  small  country  of  only  28  787  km2  surface  area  and 
around  4  500  000  inhabitants,  situated  in  the  southwest  part  of 
the  Balkan  Peninsula.  Like  the  other  Balkan  countries,  Albania 
is  located  next  to  the  subduction  boundary  between  the  African 
plate  and  the  Euro-Asiatic  plate.  This  setting  makes  the  pres¬ 
ence  of  geothermal  resources  possible.  Surface  manifestations 
of  geothermal  resources  are  found  throughout  Albania,  ranging 
from  the  region  of  Peshkopia  in  the  northeast,  where  hot  springs 
with  water  temperature  of  about  43  °C  and  inflow  above  141/s 
are  found,  through  the  central  part  of  the  country  with  different 
sources  (including  the  springs  of  Llixha-Elbasan)  with  temper¬ 
atures  above  66  °C,  to  the  Peri-Adriatic  depression  (see  Fig.  1) 
with  a  number  of  wells  (drilled  for  oil  and  gas  research)  pro¬ 
ducing  water  with  temperatures  around  40  °C  at  variable  yields. 
The  thermal  waters  in  Albania  are  only  used  for  balneology. 
This  form  of  use  dates  back  to  early  in  history,  or  to  the  time 
of  the  Roman  Empire  (i.e.  the  Sarandaporo’s  thermal  baths).  So 
far,  the  geothermal  resources  have  not  been  utilized  for  other 
purposes,  such  as  space  heating.  Estimated  temperature  mea¬ 
surements  based  on  different  geothermometers  indicates  that  the 
temperature  of  the  waters  in  the  formation  of  the  Llixha  reser¬ 
voir  may  be  above  220  °C.  The  reservoir  is  believed  to  be  in  the 
depth  interval  of  4500-5000  m.  Tirana’s  geothermal  heat  sources 
are: 

•  Underground  waters  of  the  Tirana  quaternary  depression; 

•  Underground  waters  of  the  tortonian’s  sandstones; 

•  Heat  of  the  peri-superficial  quaternary  or  tortonian  layers; 

•  Limestone’s  and  dolomites  saturated  with  artesian  waters. 

The  following  is  addressed  in  the  material: 

•  The  general  geological  conditions  in  the  areas; 

•  A  review  of  the  theoretical  basis  of  heat  transfer; 

•  Finite-volume  modeling  of  the  whole  geothermal  system  down 
to  5000  m  depth,  incorporating  both  thermal  convection  and  con¬ 
duction,  based  on  a  simple  boundary  conceptual  model; 

•  House  heating  calculations; 

•  Borehole  heat  exchangers  calculations; 

•  Economical  analyses. 

All  this  aims  at  demonstrating  that  the  thermal  water  flow¬ 
ing  from  the  Llixha  springs  is  usable  for  direct  utilization  and  the 
borehole  heat  exchanger  heating  systems  are  economically  feasi¬ 
ble,  despite  their  elevated  preliminary  cost.  This  utilization  would 
mitigate  the  electricity  supply  to  the  region,  help  improve  living 
conditions  for  the  local  community  and  protect  the  environment 
of  the  region. 

2.  Geological  background  of  the  Llixha-Elbasan  and  Tirana 
region 

2.2.  The  geological  structure  of  Llixha  region 

The  Llixha  region  is  situated  southwest  of  Elbasan.  The  region 
is  well  known  for  its  thermal  springs,  appreciated  since  ancient 
times  for  their  curative  properties.  A  geological  map  of  the 
region  and  a  representative  cross-section  (I-I)  are  presented  in 
Figs.  2  and  3. 

The  region  under  study  is  south  of  Shkumbini  river  valley. 
The  surface  relief  increases  rapidly  up  to  intermediate  eleva¬ 
tion  (300-500  m).  In  the  western  part  of  the  region,  a  system  of 
hills  declines  gradually  in  the  Cerriku  field.  The  region  has  a  rich 


hydro  system  of  small  streams  and  many  underground  water- 
systems.  The  flow  rate  for  the  underground  waters  varies  from 
100  to  200 1/h  in  Thane  up  to  5000 1/h  in  Tregan.  Generally,  the 
formations  are  composed  of  flysch  with  a  diverse  and  chaotic  mor¬ 
phology.  This  is  an  inhabited  area,  with  small  villages  clustered 
around  the  thermal  waters  no  more  than  2-3  km  from  each  other 
[1].  The  region  lies  between  two  tectonic  regions;  the  transver¬ 
sal  Vlore-Elbasan-Diber  and  the  longitudinal  Leskovik-Drini  river 
bay.  Both  of  these  connect  the  lower  part  in  connection  with 
the  northern  and  western  part  of  the  country.  In  the  context  of 
Albanian  tectonics,  the  region  represents  the  western  part  of  the 
Kruja  tectonic  zone.  The  formation  represents  the  Llixha  synclinal 
structure  limited  by  the  anticline  structures  of  Elbasan-Valeshit 
in  the  east  and  Papri  in  the  west.  The  Elbasan-Valeshi  structure 
rises  at  the  surface  through  the  calcareous  formation  of  Creta- 
Paleogenit,  while  in  Papri  the  calcareous  formations  are  at  greater 
depth.  The  orientation  of  these  structures  is  SW-NE  as  is  the  rest 
of  the  Albanides  structure  [2].  The  formations  in  the  region  are 
mainly  composed  of  flysch  (three  Oligocene  and  lower  Miocene 
(less  molasses)  sections)  while  calcareous  formations  with  lim¬ 
ited  surface  area  are  mainly  composed  of  olystolit,  through  the 
over  lined  terygen  rocks.  The  biggest  part  is  composed  of  lower 
Oligocene  (Pgs1 )  rocks,  consisting  of  6  lithological  units  as  shown 
in  Fig.  2  [1  ].  The  lowest  part  of  the  formation  is  composed  of  thin¬ 
ner  flysch  units  intercalated  by  clays  and  sandstone,  while  the 
upper  part  is  composed  of  thicker  flysch  with  conglomerates.  The 
lithology  of  the  middle  and  upper  Oligocene  formation  is  mainly 
flysch  with  differently  shaped  packs.  Interesting  is  the  lithologi¬ 
cal  evidence  of  the  upper  packs  of  the  middle  Oligocene  where 
among  the  flysches  are  located  with  lengths  up  to  100  m.  They 
are  organized  in  a  chain,  preserving  the  general  orientation  of  the 
other  strata.  Together  with  the  conglomerate  sandstone  formation 
they  form  a  chain  alternating  in  relief.  The  lithology  of  the  lower 
Miocene  (Acquitan-Burdigalian)  deposits  is  characterized  by  Mar- 
goles  alternating  with  the  clays-sandstones  as  well  as  the  presence 
of  sandstone-conglomerates. 

The  Llixha  syncline  represents  a  depressed  structure  with 
eastern  asymmetry  filled  in  the  central  part  with  terygen,  fly¬ 
sch  and  molasses  deposits.  The  eastern  part  is  distinguished 
by  an  easterly  drop.  The  tectonics  puts  the  upper  deposits 
in  contact  with  the  lower  flysch  Oligocene  formation  and  has 
made  the  surface  intrusion  of  the  Eastern  anticline  calcare¬ 
ous  formations  possible.  This  tectonic  layout  is  regional  and 
includes  the  Western  anticline  chain  of  the  Kruja  tectonic  zone 
[3]. 

The  Llixha  system  comprises  a  reservoir  which  feeds  the 
southern  part  of  the  Shkumbini  River.  The  main  hydrological  char¬ 
acteristic  of  the  region  is  the  presence  of  several  hot  springs. 
Their  position  is  related  to  the  Kruja  geothermal  zone  and 
they  are  connected  with  the  calcareous  olystolits  of  the  upper 
Oligocene  conglomerate-sandstones.  The  temperature  of  the  hot 
springs  varies,  and  ranges  from  50  to  68  °C  while  flow  rates  vary 
from  one  spring  to  another,  without  any  seasonal  characteris¬ 
tics. 


2.2.  The  geological  structure  of  Tirana 

Tirana  aquifer  is  related  with  the  syncline  deposits  of  Tirana, 
whose  morphology  represents  a  depression  of  10-12  km  of  width 
and  70-80  km  of  length,  Fig.  4  (Geological  map  of  Albania,  1984). 
The  aquifer  is  composed  of  flysch,  gravels,  sandstones,  alevrites  and 
quaternary  clays.  Well  yield  in  this  region  varies  between  7  and 
101/s,  while  for  the  Tortonians  molasses  which  are  composed  of 
sandstone,  clays  and  alevrites  the  thickness  is  60-150  m  while  the 
well  yield  is  about  3-4 1/s. 
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Fig.  1.  Map  of  Albania. 


3.  Heat  transfer  theory 

3.1.  The  differential  equation  for  heat  transfer 

This  equation  is  a  mathematical  expression  of  the  first  law  of 
thermodynamics,  the  energy  preservation  low.  The  heat  increase 


of  an  elementary  volume  AV  is  equal  to  the  thermal  energy  which 
crosses  the  surface  S.  A  solid  medium  is  considered,  which  is  not 
generating  any  energy  (so  the  energy  is  only  flowing  through  a 
surface  S).  The  temperature  T  at  a  point  P(x,  y,  z)  will  be  a  con¬ 
tinuous  function  of  the  position  and  time.  For  a  homogenous  solid 
medium,  in  which  the  thermal  volume  heat  capacity  is  independent 
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Fig.  2.  Geological  map  of  the  Llixha  region. 


of  temperature,  the  equation  is  [4]: 

2-r,  3  2T  d2T  d2T  1  dT 

dx2  dy2  dz2  a  dt 

with  a  thermal  diffusivity  of  the  solid  medium. 

In  polar  coordinates:  x  =  rcosor,  y  =  rsina,  hence: 

2j  _  3 2T  1  9 T  1  92r  92r  _  1  9T 

9r2  +  r  dr  +  r2  <)02  +  dz 2  a  dt 


Let  us  assume  that  in  point  P(x,  y,  z)  the  energy  input  per  unit  of 
time  and  unit  of  volume  is  A(x ,  y ,  z,  t),  then  Eq.  (2)  is  rewritten  as: 

(i)  V2r  +  ^3^0  =  i  3r  (3) 

X  a  dt 

If  thermal  conductivity  depends  on  position  and  temperature  the 
equation  can  be  written  as: 


d_  f  df\  d_ 
dx  y  dx  J  +  dy  y  dy  J  +  dz 


+  A 


(4) 
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Fig.  3.  Geological  cross-section  through  the  area  (I-I  in  Fig.  2). 


Its  solutions  are  often  relatively  simple  and  can  be  found  for  dif¬ 
ferent  functional  forms  of  A(x,  y,  z).  If  the  thermal  properties  are 
temperature  dependent  the  equation  is  non-linear  and  numerical 
methods  can  be  used  for  its  solution. 


3.2.2.  The  basics  hypothesis 

Let  us  describe  mathematically  the  heat  exchange  process 
between  flowing  fluids  in  a  wellbore  with  appropriate  boundary 
formations.  Let  us  assume: 
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Fig.  4.  Geological  map  of  Tirana  aquifer. 
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•  The  formation  is  homogenous  and  isotropic; 

•  The  flow  is  in  axial  direction  only; 

•  The  heat  flow  is  only  by  conduction; 

•  There  is  no  radial  temperature  gradient  in  the  wellbore; 

•  In  the  boundary  zone  of  the  wellbore  the  formation  temperature 
is  a  known  function  of  depth; 

•  The  thermal  properties  of  the  fluid  and  formation  are  constant; 

•  The  radial  temperature  in  wellbore  is  constant; 

•  The  yield  from  the  wells  is  constant; 

•  The  flow  in  wellbore  is  vertical  (ID); 

•  The  vertical  heat  conduction  -  is  much  less  than  the  horizontal 
one,  i.e.  negligible. 


conductivity  X  =  X(x,  y),  the  boundary  conditions  are: 


T(zux)  =  T\ ;  T(z2,x)  =  T2; 


-/<*</; 


dT(z,  x) 
dx 


3  T(z,x) 

=  ~rll  =  0 


(10) 


In  Table  1,  a  summary  of  the  analogies  between  hydrodynamic, 
heat  transfer,  electrostatic  and  electricity  transmissions  fields  is 
presented.  As  we  can  see  from  Table  1,  the  temperature  change 
AT  is  analogous  with  the  potential  difference  A U.  Assuming  the 
length  unit,  lSl  the  dimensionless  coordinates  are: 


3.2.  Initial  and  boundary  conditions 

In  order  to  solve  the  thermal  diffusion  equation  let  us  assume 
that  the  formation  temperature  is  function  of  the  position  (x,  y,  z) 
and  time  (t)  [4].  In  addition  initial  and  boundary  conditions  must 
be  specified.  A  specific  moment  of  time  is  chosen  as  the  origin  of 
the  time  coordinate.  At  that  time  the  temperature  distribution  is: 

T(x,y,  z,t  =  0)  =  f(x,y,  z)  or  T(r,&,z,t  =  0)=f(r,0,z)  (5) 


The  relationship  X  =  X(x,  z)  is  analogous  with  p  =  p{x,  z).  A  partic¬ 
ular  care  should  be  used  in  unstable  field  modeling.  Considering 
the  hydrodynamics  analogy  with  the  thermal  one,  in  particular 
the  transitory  flows  of  an  incompressible  fluid  through  a  porous 
medium,  the  main  characteristics  of  this  regime  are  the  high  values 
of  the  hydraulic  diffusivity  [2]: 


If  the  radials  symmetric  case  of  a  flowing  wellbore  is  considered, 
as  an  example.  After  a  certain  production  time  the  well  is  shut 
down.  To  determine  the  temperature  distribution  in  the  wellbore 
during  the  shutdown  time,  the  end  of  production  is  considered  as 
the  origin  of  the  time  coordinate.  In  this  case  the  temperature  /(r, 
z)  must  be  known,  i.e.  the  initial  conditions.  To  specify  the  tem¬ 
perature  field  of  a  medium  the  boundary  conditions  must  also  be 
known  beforehand.  Different  kinds  of  boundary  conditions  can  be 
specified: 


(12) 


Its  thermal  field  analogue  is  the  thermal  diffusion  coefficient  (a). 
The  analogies  principle  can  only  be  used  in  the  case  when  the 
dimensionless  times  are  equal.  In  thermal  field  modeling,  the  fact 
that  the  temperature  is  not  affected  by  surface  topography  and 
groundwater  flow  is  very  important.  To  determine  the  boundary 
conditions  temperature  maps  at  depth  are  widely  used  [5]. 


•  Surface  temperature  is  known.  It  can  be  constant  or  function  of 

position  and  time,  T =/(x,  y,  z,  t). 

•  The  amount  of  energy  flowing  through  the  surface  is  known: 

dT  f  d  1 

qs(x, y,  z,  t)  =  -X  —  <  —  =  the  derivation  perpendicular  to  the  surface  >  (6) 

dn  [  dn  J 

•  Linear  surface  heat  flow.  In  such  a  case  the  amount  of  energy 
transmitted  through  a  given  surface  is  proportional  with  the  tem¬ 
perature  difference  between  the  surfaces  and  the  surroundings. 


3.4.  The  unstable  temperature  field 

Before  starting  investment  in  a  geothermal  project,  the  stability 
of  the  temperature  field  involved,  in  space  and  time,  needs  to  be 
confirmed  as  well  as  the  projects  overall  sustainability.  The  meth¬ 
ods  that  can  be  used  to  answer  such  questions  are  numerous,  but 
here  we  will  apply  the  finite  element  method  [6].  Considering  the 
differential  equation: 


qs=a(Ts-T0)  (7) 

where  Ts  is  the  surface  temperature;  T0  is  the  surrounding  tem¬ 
perature;  a  is  the  heat  transfer  coefficient,  and 

lim^ooTs  =  T0  (8) 


d_ 

dx 


+  q(x,y,  z)  — 


=  0 
(13) 


The  functional  of  this  equation  is  [6]: 


•  The  connecting  surface  between  two  media  with  conductivities 
A-t  and  X2,  respectively.  If  T ^  and  T2  are  the  temperatures  of  the  X  =  Xy  +Xr  = 
media,  then 


kx  ( dT\2  ky  ( dT 
2  ydxj  2  dy 


2 


Tj  Is  =  T2|s 


dTx 

dn 


s 


x  dxdydz + 


dxdydz 


(14) 


3.3.  Temperature  field  modeling 


Let  us  divide  this  functional  into  two  parts: 


In  order  to  create  a  physical  model  to  enable  the  stable  state  of 
heat  flow  through  a  formation  to  be  studied,  the  analogy  princi¬ 
ple  may  be  used.  Let  us  assume  that  n  dimensionless  parameters 
/?,. .  .dn  describe  the  heat  conditions  in  the  formation  and  sur¬ 
roundings.  In  such  a  case,  the  model  results  should  be  expressed 
as  n  dimensionless  parameters  /j , . . . ,  Tn.  The  conditions  that  must 
be  fulfilled  in  order  to  use  the  model  in  a  real  reservoir  are  l\  = 
/j , . . . ,  / 1  =  Tn.  For  stable  heat  transmission,  the  analogy  with  the 
electrostatic  field  can  be  used.  For  the  2D  case  and  with  the  thermal 


To  integrate  the  Eq.  (14)  let  us  assume  the  time  interval  (t,  t+  At) 
through  into  dTjdt  =  C 
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Table  1 


The  analogy  between  different  physical  fields  (based  on  [4]). 

Hydrodynamics 

Heat  transfer 

Electrostatics 

Electricity  transmission 

Pressure 

Temperature 

Electrostatic  potential 

Potential 

P 

T 

0 

U 

Pressure  gradient 

Temperature  gradient 

Field  vector 

Potential  gradient 

-vP 

-VT 

E=-V0 

-VU 

Permeability/viscosity  /</|x 

Thermal  conductivity  X 

Dielectric  constant  e\4n 

Specific  resistance  p 

Velocity  vector  v  —  --jjVp 

Heat  amount  transfer  q  =  -  VT 

Dielectric  displ.  (e/47r)lf  =  -(e/47r)V0 

Current  /  =  -  VU 

Isobaric  surface  p  =  C 

Isothermal  surface  T=C 

Surface  with  electrostatic  potential  0  =  C 

Surface  with  potential  s  =  C 

Impermeable  layer  =  0 

Split  surface  ^  =  0 

Force  line  ^  =  0 

Potential  line  =  0 

7):  temperature  values  in  the  nodes  7=1,4 


ax 


dTi  dTi 


dXv  dXj 
+ 


dTi  ’ 


T  = 


6A 


{ Ni  N2  N3 


(  T,  ' 

N4 }  <  l2 


h 

[t4) 


where 


(16) 


N/(x,  y,  z)  =  (a,  +  b{x  +  qy  +  d,z);  7  =  1,4 

The  partial  derivates  of  the  functional  can  be  calculated  as  follows: 

ax 


4.  District  heating  system  using  the  geothermal  water 

In  general  the  geothermal  water  for  district  heating  systems  is 
taken  directly  from  low  temperature  reservoirs.  Another  way  is  to 
use  the  geothermal  water  through  the  heat  exchangers  to  heat  up 
the  fresh  water.  The  hot  water  can  be  stored  in  tanks  if  appropriate. 
This  water  then  is  transmitted  to  the  buildings  and  can  be  used  for 
heating  and  as  tap  water.  Heat  flow  to  the  buildings  is  controlled  by 


ax 

Wi 


=  < 


art 

ax 

W2 

ax 

ar3 

ax 


>  = 


l  dT 


4  2 


/ 

b] 

b\b2 

b\b2 

b\b4 

c2 

L1 

ctc2 

ctc3 

C1C4 

di 

d\d2 

d\  d3 

d\d4 

\ 

h 

kx 

b2b\ 

b2 

u2 

b2b3 

b2b4 

ky 

+ 

36  A 

c2ct 

c2 

l2 

C2C3 

C2C4 

+  kz 

36  A 

d2d\ 

d\ 

d2d3 

d2d4 

1 

T2 

36  A 

b2b\ 

b3b2 

bl 

b2b4 

C3Ct 

C3C2 

C2 

l3 

C3C4 

d3d\ 

d3d2 

d2 

d3d4 

1 

t3 

{ 

Mi 

b/[b2 

b4b3 

b 2 
u4 

C4C1 

C4C2 

C4C3 

C2 

l4 

d4d\ 

d4d2 

d4d3 

d2 

u4 

J 

t4 

/  /  /  N?dxdydz 

J  J  Jn 

/  /  /  NiN2dxdydz 
JJJn 

/  /  /  NiN3dxdydz 
JJJn 

/  /  /  NiN4dxdydz 
JJJn 

/  /  /  N2N-idxdvdz 
JJJn 

/  /  /  Nldxdydz 

JJJn 

/  /  /  N2N3dxdydz 
JJJn 

/  /  /  N2N4dxdydz 
JJJn 

/  /  /  N3N^dxdydz 

J  J  J  n 

/  /  /  N3N2dxdydz 
JJJn 

/  /  /  Nldxdydz 

J  J  Jn  3 

/  /  /  N3N4dxdydz 
JJJn 

/  /  /  N4Nidxdydz 

J  J  J  n 

/  /  /  N4N2dxdydz 

J  J  J  n 

/  /  /  N4N3dxdydz 

J  J  J  n 

/  /  /  Njdxdydz 

JJJn 

dT2 

1 

0 

dt 

Q.A 

1 

0 

dT3 

4 

1 

0 

dt 

1 

0 

(17) 


Let  us  make  the  following  replacements: 


kx 

b\  M2  M3  M4 

Mi  b\  M3  b2b4 

,  ky 

36A 

C2  C\C2  C\C3  C\C4 

C2C\  c2c3  C2C4 

+  kz 

36A 

d\  did2  M3  d,d4 

d2d\  d 2  d2d3  d2d4 

36  A 

b3bx  b3b2  b2  b3b4 

c3ct  c3c2  c2  c3c4 

d3d\  d3d2  d2  d3d4 

Mi  M2  M3  b\ 

C4C\  C4C2  C4C3  c2 

d4d\  d4d2  d4d3  d^ 

=  [H] 


(18) 


/  /  /  N?dxdydz 
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/  /  /  N-[N2dxdydz 

J  J  J  Q 
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JJJn 
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JJJn 

/  /  /  N2N^dxdydz 

J  J  J  n 

/  /  /  Nldxdydz 

JJJn 
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JJJn 
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JJJn 

/  /  /  N3N2dxdydz 

J  J  J  P2 

/  /  /  Nldxdydz 

JJJn 

/  /  /  N3N4dxdydz 
JJJn 

/  /  /  N4N^dxdydz 

J  J  J 

/  /  /  N4N2dxdydz 
JJJn 

/  /  /  N4N3dxdydz 
JJJn 

/  /  /  Njdxdydz 

JJJn 

=  [p] 


(19) 


QA 


=  [F] 


(20) 


where  [H]  is  the  matrix  of  conductivity  of  the  thermal  field,  [P]  is 
the  matrix  of  instability  of  the  thermal  field,  [F]  is  the  source  vector. 

After  these  replacements  the  definitive  form  of  Eq.  (17) 
becomes: 


[H] 


+  [p'|-|fW 


=  0 


(21) 


the  mass  flow.  In  the  following  are  given  the  basics  calculations  for 
the  geothermal  district  heating.  The  main  elements  of  the  geother¬ 
mal  district  heating  are  the  radiators  and  also  this  use  is  affected 
by  the  water  thermal  energy,  building  heat  loss,  pipe  heat  loss  and 
building  energy  storage.  A  short  description  for  all  of  them  is  given 
below. 
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Table  2 

Estimated  water  temperature  at  depth  based  on  different  geo-thermometers. 


Geothermometer 

Spring  “Nosi”  Llixha,  Elbasan  (°C) 

Fournier 

254 

Truesdell 

235 

Na  +  Ka  +  Ca 

143 

20' S' 


* 


A 


Fig.  5.  Geomorphologic  map  of  the  Llixha  region. 


4.2.  Calculation  procedure  for  house  heating  systems 

Radiators  are  the  heat  exchangers  and  they  make  the  trans¬ 
fer  of  the  heat  from  the  geothermal  water  to  the  indoor  possible. 


•  bottom_cent£ 

•  bottom'surro 
top_center 

•  top"surround 


Stat  ic 
T  emperature 
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Fig.  7.  The  temperature  magnitude. 
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Position  (m) 

Fig.  8.  The  density  magnitude. 


The  relative  heat  capacity  of  the  radiator  is  given  by  the  following 
relation: 


Qrad  =  (  ATm  \4/3 

%ad  VATm(J 


(22) 


where  the  subscript  0  denotes  the  design  conditions.  The  temper¬ 
ature  difference  A Tm  is  calculated  [7]: 


(Ts  -Tj)  —  (7V  —  Tj)  =  Ts-Tr 

In  ((Ts  -  T,)/(Jr  -  Tj))  In  ((Ts  -  T,)/(Tr  -  Tj)) 


(23) 


The  water  thermal  energy  which  the  geothermal  water  will  give  to 
the  indoors,  through  the  radiators  is: 

Qrad  =  cpm(Ts  -  Tr)  (24) 


The  relative  thermal  energy  given  by  the  geothermal  water  is  [7]: 

Qrad  =  m(Ts  -  Tr) 

Q^ad  mo(Tso  -  T,o) 


•  botl  om^cente 

•  botlom'surro 
top_center 

•  totTsurround 


Velocity 
Nagnit  ude 
(m/s) 


Fig.  6.  The  finite  volume  grid. 


Fig.  9.  The  velocity  magnitude. 
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Table  3 

Thermal  characteristics  of  the  rocks. 


Age 

Lithology 

Thermal  conductivity  (W/m°C) 

Thermal  resistance  (m°C/W) 

Tortonian 

Sandstone 

Clay 

1.6 

1.5 

0.625 

0.666 

Table  4 

Physical  characteristic  of  the  sandstones. 

Layer 

Filtration  coefficient  Water  conductivity  (m/day) 

Specific  yields  (lm/s) 

Tortonian  sandstone,  upper  part 
Tortonian  sandstone,  lower  part 

0.08  3.2 

0.08  5.3 

0.038 

0.038 

The  biggest  part  of  the  energy  loss  in  such  system  may  come  from 
the  buildings.  These  losses  can  be  calculated  as  [7]: 

Qioss  =  ki(Ji  -  To)  (26) 

where  building  heat  loss  factor  /</  is  constant.  The  relative  losses  are 
[7]: 

Qioss  =  Tj  -  Tp 
QlossO  Tjq  -  Too 

The  other  losses  in  such  system  are  the  pipe  heat  losses.  To  deter¬ 
mine  the  amount  of  losses  in  the  pipes  it  is  necessary  to  know  the 
pipe  transmission  effectiveness  parameter  r.  Its  value  [17]  is  given 
by: 


A  building’s  heat  storage  (dependent  on  the  amount  of  the  thermo¬ 
insulation)  can  be  very  helpful  for  the  heating  system.  Normally  the 
builds  are  like  a  heat  store.  The  building  heat  storage  is  calculated 
by  [7]: 


w  ~  zQnet  ~  r(Qsup 


Qioss )  —  —  fnCp(Ts  —  Tr)  —  k/(Tj  —  To) 


(33) 


5.  Borehole  heat  exchangers 

5.2.  Equation  for  heat  flow 


x  = 


Ts 

Ti 


Jg 

Tg 


=  exp 


(28) 


The  reference  value  of  r  can  be  calculated  from  the  reference  flow 
conditions  as  follows: 


T0 


Tso  -  Tgo 
Tio  -  Tgo 


=  exp 


Up 

m0cp 


(29) 


The  values  of  Up  and  cp  are  assumed  to  be  constant  in  the  system. 
By  combining  Eqs.  (28)  and  (29),  the  transmission  effectiveness  can 
be  calculated: 


r  =  x™o/m  (30) 

The  supply  temperature  to  the  house  can  be  calculated: 

Ts  =  Tg  +  (T,  -  Tg)r  =  Tg  +  (T,  -  Tg)x™°lm  (31) 

and  the  return  water  temperature  at  the  pumping  station  is  calcu¬ 
lated: 


Ti  =  Tg  +  (Tr  -  Tg)r  =  Tg  +  (Tr  -  Tg)r™°lm 


(32) 


8.51  e-  D8 
8.  D9e-  08 
7.66e-08 
7.246-08 
6.81  e-  08 
6.39e-  08 
5.96e-D8 
5.53e-  08 
5.1  le-08 
4.68e-  08 
4.26e-  08 
3.83e-  08 
3.41  e-  08 
2.98e-  08 
2.55e-  08 
2. 1 3e-  08 
1 .70e- 08 
1 .28e- 08 
8.52e-  09 
4.26e-  09 
4.53e-  1  2 


The  equation  for  rate  of  heat  transfer  from  the  fluid  in  the  heat 
exchanger  to  the  earth  mass  is: 

^  =  Ux  AT  (34) 

where  Q.  is  the  rate  of  heat  exchanger  (BTU/h  or  W)  for  the  whole 
heat  exchanger  length;  L  is  the  length  of  heat  exchanger  (m);  U  is 
the  conductance  rate  for  heat  transfer  from  the  circulating  fluid  to 
the  earth  (BTU/h/°F  or  W/°C/m)  for  the  operating  conditions:  AT  is 
the  difference  in  fluid  temperature: 

AT  =- 2-2E  -  t0  (35) 

T0  is  the  earth  temperature  (°F  or°C);  T\  is  the  fluid  entry  temper¬ 
ature  (°F  or°C);  T2  is  the  fluid  exit  temperature  (°F  or°C). 


Fig.  10.  The  velocity  vectors. 


Fig.  11.  Calculations  for  the  radiators. 
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Lithological  column  of  the  well  Nr. 
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Fig.  12.  Lithological  column  of  a  well  drilled  in  Tirana. 


The  conductance  term  U  for  heat  flow  from  fluid  in  the  heat 
exchanger  to  the  earth  can  be  estimated  with  the  conductance 
coefficient  for  composite  cylinders.  The  impedance  to  heat  flow 
is  caused  by  the  thermal  resistance  of  the  pipe  wall  and  the  soil 
cylinder  around  the  casing.  Fluid  surface  resistance  films  are  small 
relative  to  the  other  terms  and  are  encompassed  in  the  two  resis¬ 
tance  terms.  We  can  express  U  as: 


U  = 


2tt 

soil  resistance  +  pipe  resistance 


2n 

Rs  +  Rp 


(36) 


show  that  the  water  of  these  hot  springs  can  be  used  for  the  district 
heating  of  the  village  community. 


6.2.  The  3-D  modeling  of  the  temperature  field  ofLlixha,  Elbasan 
region 

A  finite  volume  model  was  set  up  for  a  crustal  volume  with  an 
area  of  10  km  x  10  km  and  5  km  thickness  to  model  the  temper¬ 
ature,  density  and  fluid  velocity  distribution  in  the  Llixha  region. 
The  grid  is  shown  in  Fig.  6.  Here  it  is  assumed  that  the  medium 
is  homogeneous  and  isotropic  and  that  kx  =  ky  =  kz  =  2  W/mK  [9]. 
We  also  know  that  Q  =  201/s  (corresponding  to  mj  =  Q/ 6  =  3.3 1/s  or 
3.224 kg/s  for  each  of  the  hot  springs),  cv  =4180J/kg°C  [8].  The 
temperature  at  depth  in  the  formation  is  set  at  221  °C  while  the 
temperature  of  the  water  at  the  surface  is  in  the  range  60-65  °C. 
The  temperature  gradient  of  the  surroundings  is  assumed  1 2  °C/km. 
The  modeling  software  FLUENT  is  used  to  solve  the  problem,  it  pro¬ 
vides  calculation  results  for  temperature,  density  and  velocity  for 
the  volume  modeled.  In  the  model  water  flows  with  a  velocity  of 
1.25  x  10~7  m/s.  The  results  for  temperature,  density  and  velocity, 
as  well  as  velocity  vectors,  are  shown  in  Figs.  7-10. 


7.  Basic  district  heating  design 

7. 2 .  Calculation  for  the  radiators 

The  steps  used  to  calculate  the  radiator  parameters  are  given 
in  Section  4.  In  these  calculations,  it  was  assumed  that  the  sup¬ 
ply  temperature  of  the  water  is  60  °C,  the  reference  temperature 
65  °C,  and  the  ground  temperature  6  °C.  The  calculations  were  done 
for  4  different  scenarios:  The  indoor  temperature  is  assumed  to 
change  at  the  range  of  18-20°C,  the  outdoor  temperature  in  the 
range  of  - 1 0  to  -4  °C,  the  return  water  temperature  in  the  range  of 
33-40  °C,  the  reference  inflow  in  the  range  of  3.93-7  kg/s,  the  ref¬ 
erence  system  inflow  18-24  kg/s.  Based  on  these  data,  the  relative 
heating  of  the  radiators  was  calculated  as  0.83-0.92,  the  relative 
heating  of  the  building  as  0.87-0.88,  and  the  transmissivity  coeffi¬ 
cient  r  =  0.94.  To  be  within  these  parameters,  it  is  sufficient  that  the 
supply  water  temperature  be  60  °C,  and  the  inflow  of  the  system 
1 6-27  kg/s.  Thus  the  parameters  for  the  Llixha  thermal  springs  sat¬ 
isfy  all  these  demands.  Fig.  1 1  shows  the  results  of  the  calculations 
for  the  radiators  supply  and  return  water  temperatures  vs.  indoor 
temperature. 


7.2.  Borehole  heat  exchanger  calculations 


6.  Results  of  simple  hot  spring  modeling 

The  six  hot  water  springs  at  Llixha  in  the  Elbasan  region,  have 
water  temperature  up  to  65  °C  and  flow  rate  up  to  231/s.  Esti¬ 
mated  reservoir  temperature  at  depth  associated  with  the  hot 
springs,  based  on  the  chemical  composition  using  different  geo¬ 
thermometers,  is  given  in  Table  2. 

These  values  show  that  the  water  is  coming  from  great  depth 
where  the  average  temperature  is  over  200  °C.  As  it  is  shown  in 
Fig.  5  [8  ],  the  hot  springs  are  situated  in  the  middle  of  the  village.  The 
hot  water  has  been  used  only  for  balneology  for  several  centuries 
possibly  since  the  time  of  the  Roman  Empire.  The  first  modern  use 
dates  back  to  1937  with  the  building  of  the  “Hotel  Park”  medical 
centre.  The  use  of  the  water  flowing  from  these  springs  can  help 
to  improve  the  economical  effectiveness  district  heating  in  the  vil¬ 
lage.  The  realization  of  such  a  project  will  allow  the  utilization  of 
geothermal  water  as  an  energy  resource  for  the  first  time  in  Alba¬ 
nia.  The  purpose  of  the  calculations  presented  in  this  chapter  is  to 


Heat  exchanger  length  is  calculated  through  Eq.  (34).  The  ther¬ 
mal  characteristics  of  the  rocks  drilled  by  the  well  of  Fig.  12  are 
given  in  Table  3,  while  the  physical  characteristics  of  Tirana  torto- 
nian  sandstones  are  given  in  Table  4. 

For  this  well  is  known: 


Dynamic  level 

100  m; 

Static  level 

50m; 

Level  falls 

40  m; 

Pump  depth 

108  m; 

Average  yield 

ll/s; 

pH 

7.5-8; 

Water  temperature 

13-15°C; 

General  strength 

10-15  °Gj; 

Mineralization 

250-350  mg/1; 

Formula 

HCOs-Ca-Mg; 

The  thermal  resistance  for  the  formations  Rp  drilled  by  the  well 
is  calculated  as  follows: 


(40  x0.666  +  34x0.625) 


0.647 


m  x°  C 
W 


40  +  35 


2544 
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Fig.  13.  The  economical  analyses  of  three  different  systems. 


The  conductance  rate  U  is  calculated: 


U  = 


2tt 

¥ 


2tt 

0.647 


9.711 


W 

mx°C 


According  to  the  measurements,  the  earth  temperature  in  100  m 
depth  in  Tirana  is  r0  =  18°C.  The  fluid  exit  temperature  is  3.5  °C 
lower  than  the  fluid  entry  temperature.  The  difference  temperature 
of  the  system  is  calculated: 


AT  =  T0  -  Tl  ty  2  =  6.75  °C 

For  these  parameters  the  heat  exchangers,  for  the  installed  capacity 
of  lOOkW  result: 


L  = 


d 

U  x  AT 


100  x 103 
9.711  x  6.75 


=  1525m 


So  the  problem  is  solved  through  drilling  process  of  1 5  wells  1 00  m 
depth  each.  The  economic  analyses  (profit-expenditure)  of  three 
systems:  Water-Water  geothermal  pumps,  Earth-Water  geother¬ 
mal  pump  and  the  existing  coal  heater  are  presented  in  Fig.  13.  So  it 
can  be  clearly  seen  that  the  maximal  profit  is  for  the  Water-Water 
system  than  for  the  Earth-Water  and  finally  for  the  coal  heater.  The 
payback  period  for  the  geothermal  heating  pumps  varies  from  2.8 
up  to  10.75  years. 


8.  Conclusions 

Based  on  the  calculations  presented  in  this  report  the  following 

can  be  concluded  regarding  utilization  of  the  hot  springs  of  the 

Llixha-Elbasan  hot-spring  area  in  Albania: 

•  The  water  temperature  is  expected  to  be  stable  in  the  future; 

•  The  geothermal  water  from  the  Llixha  hot  springs  fulfils  all 
requirements  for  district  heating’s  in  the  region; 

•  Considerably  higher  temperature  expected  through  further  well 
drilling; 

•  Albanian  geothermal  regime  allows  different  scale  borehole  heat 
exchangers  applications; 

•  Demographic  and  geological  features  of  the  student’s  city  allow, 
and  furthermore  are  feasible  the  borehole  heat  exchanger’s  uti¬ 
lization. 
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